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ABSTRACT 



Context. Coronagraphy is a powerful technique to achieve high contrast imaging and hence to image faint companions around bright 

targets. Various concepts have been used in the visible and near-infrared regimes, while coronagraphic applications in the mid-infrared 

remain nowadays largely unexplored. Vector vortex phase masks based on concentric subwavelength gratings show great promise for 

such applications. 

Aims. We aim at producing and validating the first high-performance broadband focal plane phase mask coronagraphs for applications 

in the mid-infrared regime, and in particular the L band with a fractional bandwidth of ~16% (3.5-4.1 /im). 

Methods. Based on rigorous coupled wave analysis, we designed an annular groove phase mask (AGPM) producing a vortex eff'ect in 

the L band, and etched it onto a series of diamond substrates. The grating parameters were measured by means of scanning electron 

microscopy. The resulting components were then tested on a mid-infrared coronagraphic test bench. 

Results. A broadband raw null depth of 2 x 10~^ was obtained for our best L-band AGPM after only a few iterations between design 

and manufacturing. This corresponds to a raw contrast of about 6 X 10~^ (10.5 mag) at 2A/D. This result is fully in line with our 

projections based on rigorous coupled wave analysis modeling, using the measured grating parameters. The sensitivity to tilt and 

focus has also been evaluated. 

Conclusions. After years of technological developments, mid-infrared vector vortex coronagraphs finally become a reality and live 

up to our expectations. Based on their measured performance, our L-band AGPMs are now ready to open a new parameter space in 

exoplanet imaging at major ground-based observatories. 

Key words. Instrumentation: high angular resolution, coronagraphy - Stars: planetary systems 



1. Introduction 

Direct imaging of exoplanets and close environments of stars has 
recently made a long-awaited breakthrough ([Marois et al.|2QQ8[ 



Kalas et al.|2Q08l [Lagrange et al.|20Q9a||2Q10[|Absil & Mawet| ( |Mawet et al 



2Q10| ). While coronagraphs have sometimes been parts of these 
discoveries, their role was restricted to merely mitigate detec- 
tor saturation. The reason for the still limited impact of stellar 
coronagraphs in direct imaging is twofold. First, the wavefront 
quality and stability provided by current instruments on most 
ground-based telescopes at short wavelengths (up to the near- 
infrared) is not high enough (Strehl ratio of the order of 50%). 
Second, the coronagraphic devices offered until a few years ago 
date back to the original century-old invention by French as- 
tronomer Bernard Lyot (Lyot 1939 ). Paradoxically, the diversity 
of new coronagraphs developed in the lab during past 15 years is 
far greater than the number of real science discoveries that they 
have enabled so far ( ,Guyon et al. 2006,). It is also fair to say 
that, apart from very few exceptions ([Riaud et al.|20Q6[|Mawet 



tests and operations ( [Mawet et al.||2012] ), from which so much 
can be learned while enabling limited but actual science despite 
the average image quality in the near-infrared. 

As a few dedicated experiments have already suggested 
201Q[ [Serabyn et"aLl|2010| ), this situation is 



about to change. Indeed, second-generation instruments start to 
roll in at major observatories, e.g., the Gemini Planet Imager 
(GPI, [Macintosh et al. 2008 ), VLT/SPHERE ( Beuzit et al.|20081 
Kasper et al.|2012 ), and a few others. They promise far greater 



let al."2009l |20TT| [Serabyn et al. 



2009b i Quanz et ai.||2010t [Boccaletti et al.| 2012), investments 



2009 2010; "Lagrange et al. 



have mainly been made to develop the technologies in the lab 
rather than actually bring them to the telescope for actual on- sky 



image quality and stability, enabling new-generation corona- 
graphs to deliver their best performances. Not to be forgotten 
are first generation instruments, which still possess untapped po- 
tential, that only 10 year s or so of operation and understanding 
allow us to fully unleash ( jGirard et al.|2012[|Mawet et al.|2012| ). 
Here we present the successful outcome of 8 years of techno- 
logical developments to mak e the so-c alled annular groove phase 
mask coronagraph (AGPM, jMawet e t al. 2005a) a reality ready 
to be installed at the telescope. The AGPM is an optical vortex 
made out of a subwavelength grating. The vortex coronagraph is 
one of the most advanced new-generation coronagraphs and its 
interest lies in its ability to reach high contrast at very small inner 
working angles (IWA), while maintaining high throughput over 
a full 360° field of view. This has already been proved in the liter- 



C. Delacroix et al.: Laboratory demonstration of a mid-infrared AGPM vector vortex coronagraph 



ature with several suc cessful demonstrations and results, b oth in 
the lab and on sky (Mav^et et al.|2010[[Serabyn et al.|2QlQ| ) but at 
shorter wavelengths (visible and near-infrared). This paper tack- 
les the difficult problem of designing and manufacturing vector 
vortex coronagraphs for the mid-infrared, where the demand for 
efficient coronagraphs is increasing following the recent success 
of high contrast imaging of exoplanets and circumstellar d isks in 
the L band around 4 jum (Moerchen et al._2QQ7i Lagrange et al. 



|2QTQl[QlIi5^ et al.|201Q||2QlH|Kenworthy et al.|2Q13| ). L band 
is indeed an ideal filter for exoplanet searches on ground-based 
telescopes: 

- The L-band contrast of planetary-mass companions with re- 
spect to their host stars is predicted to be more favorable than 
at shorter wavelengths (Baraffe et al. 2003} Fortney et al. 
2008 [ [Spiegel & Burrows 2012| ) so "that lower-mass, older 



objects can be addressed. 
- An additional advantage of L band over shorter wavelengths 
is the better adaptive optics image quality, with Strehl ratios 
typically between 70% and 90%. 

These advantages by far compensate the increased sky back- 
ground in the thermal infrared. In this context, the AGPM will 
be competitive, for it is designed to reduce the stellar contri- 
bution by AL' > 7.5 magnitudes at very small IWA (down to 
0.9A/D). An L-band AGPM on an AO-assisted 10-m class tele- 
scope should be capable of imaging giant planets at a projected 
separation of only about 1 AU from stars located at 1 Ope. 

After reviewing the principle of vortex coronagraphs and of 
the AGPM in ^ we describe in ^ the design, manufactur- 
ing and expected performance of L-band AGPM coronagraphs. 
In ^ we detail the optical setup of the coronagraphic test- 
bench and present the measured performance of our best L-band 
AGPM. Finally, we conclude by discussing the perspectives of 
on- sky applications and of future improvements. 



2. The AGPM: a vector vortex coronagraph 

Optical vortices occur when the phase structure of light is af- 
fected by a helical ramp around the optical axis, e^^^, where 6 
is the focal plane azimuthal coordinate and / is the vortex topo- 
logical charge, i.e., the number of times the geometric phase (pp 
( |Pancharatnam|| 1 95^ see also below) accumulates 27i along a 
closed path s surrounding the phase singularity: 



/: 



— &)^^pds 



(1) 



The phase dislocation forces the amplitude to zero at its cen- 
tre, which is a singularity. Nature indeed prevents the phase 
from having an infinite number of values at a single point in 
space, which is non-physical, by simply nulling the light locally. 
When centred on the diff'raction pattern of a star seen through a 
telescope, optical vortices aff'ect the subsequent propagation by 
redirecting the on-axis starlight outside the geometric image of 
the pupil. This diff'racted starlight is then blocked when pass- 
ing through a diaphragm in the so-called Lyot stop (LS) plane, 
slightly undersized compared to the input pupil. 

The perfect starlight attenuation of an optical vortex coron- 
agraph was proven analytically (Mawe ret al.||2005a[ Foo et al. 
,2005 i Jenkins 2008 ) for any non-zero even values of / by apply- 
ing the result of the Webe r- Schafheitlin integraQ ( [Abramowitz 
& Stegun|1972 Sneddon||1951J to the Fourier transform of the 




Fig. 1. Scanning electron microscope (SEM) picture of the cen- 
tre of an Annular Groove Phase Mask (AGPM) made out of dia- 
mond and dedicated to coronagraphic applications in the L band. 



product of a vortex ramp phase {e^^^) with an ideal Airy pattern 
( ^y^ ), where R is the circular input pupil radius, k the wave 
number, and p the radial coordinate in the focal plane. For in- 
stance, with a topological charge / = 2, the amplitude of the 
electric field in the Lyot stop plane, with (r, y) coordinates, be- 
comes 



^Ls(^r) 



, r<R 



(2) 



This integral reduces to Sonine's for the / = 2 case. 



The signal is perfectly nulled inside the pupil area, while the 
redirected light outside the pupil area is blocked by the Lyot stop. 
Naturally, this occurs only for the on-axis starlight whereas light 
from off'- axis objects misses the centre of the vortex and prop- 
agates normally. The vortex coronagraph is known for its near- 
perfect coronagraphic properties: small IWA (down to 0.9/1/Z)), 
high throughput, clear off'-axis 360° discovery space, and sim- 
plicity. 

An optical vortex can be either scalar or vectorial. Scalar 
optical vortices are based on a structural phase ramp such as 
a piece of glass shaped as a helix. The accuracy to be achieved 
in the glass thickness is a real technological challenge which 
still remains poorly mastered, even using the most recent micro- 
fabrication techniques. More importantly, the induced phase pro- 
file of scalar vortices is highly dependent on the wavelength of 
the incoming light. Therefore, these are not suitable for sig- 
nificantly wide spectral bands such as those commonly used 
as astronomical filters. A vectorial (or vector) vortex corona- 
graph (VVC), on the other hand, can overcome these limitations. 
Instead of a structural phase ramp, the VVC is based on a circu- 
larly symmetric halfwave plate, i.e., a mask that aff'ects the trans- 
verse polarisation state cyclicly around the centre. This creates 
a geometric phase (pp known as the Pancharatnam-Berry phase, 
which is half of the solid angle subtended by the polarisation cy- 
cle on the Poincare sphere (Pancharatnam 1956; Berry 1987). It 
was rigorously shown (Mawet et al. 2005a ; Niv et al. 2006) that 
vectorial vortices present the same phase ramp as scalar vortices. 

Diff'erent types of VVCs exist. A technology using liquid 
crystal p olymers (LCP-VVC) recently sh owed excellent results 
(Serabyii eFaLllMTQl [Mawet et aL][2QTQl ). Yet, this technology 
is limited so far to the visible and near-infrared wavelength 
regimes (up to 2.4 /im), and is not suitable for the highly cov- 
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eted mid-infrared region. Another technological route is the 
Annular Groove Phase Mask (AGPM, see Fig. [T]) using sub- 
wavelength gratings (S G-VVC), which are pa r ticularly adapted 
to l onger wavelengths ( |Mawet et al. 2QQ5b[ [Delacroix et al. 
|2Q1 Q). Sub wavelength gratings (SGs) are gratings with a period 
A smaller than the illuminating wavelength A divided by the re- 
fractive index n of the substrate (assuming that the incident light 
is perpendicular to the grating surface and that the surrounding 
medium is air). Thanks to their property of form birefringence, 
which gives two different refractive indices hje (transverse elec- 
tric) and ^TM (transverse magnetic), such SGs can produce a 
phase shift between the two polarisation components: 



AOte-tm(^) 



27r, 



-/zA^form(/^) 



(3) 



By carefully selecting the grating parameters (substrate and ge- 
ometry), this phase shift can be made quasi-independent of the 
wavelength over a wide spectral bandwidth, there by synthesiz- 
ing an artificial birefringent achromatic wave plat e ([Kikuta et al. 
|1997| ). We have recently shown (Delacroix et al. 2012b) that dia- 
mond is a good material for manufacturing achromatic half-wave 
plates for mid-infrared wavelengths (e.g., L and N bands centred 
respectively around 3.8 yum and 10 jim). Here, we focus on ap- 
plications in the L band, ranging from 3.5 to 4.1 /urn. 

3. AGPM-L design and manufacturing 

3.1. Design 

The design of the grating was conducted in com plete synergy 
with the manufacturing process, described in Sect. 32_ We per- 
formed realistic numerical simulations using the rigorous cou- 
pled wave analysis (RCWA), which resolves the Maxwell equa- 
tions in the frequency domain and gives the entire diffractive 
characteristics of the studied structure (Mohar am & Gaylord| 
|1981| . Although the theory states that a perfectly achromatic 
coronagraph would provide a perfect nulling of the on-axis 
starlight, in practice, imperfections remain. In this paper, we use 
the null depth metrics instead of the peak-to-peak attenuation. 
The null depth quantifies the contrast provided by the corona- 
graph, integrated over the whole point spread function (PSF). As 
explained further in Sect. |4.3| this metrics is well suited to coro- 
nagraphy since it takes into account the changes induced by the 
coronagraph itself to the PSF profile. The theoretical null depth 
A^theo is given by the f ollowing wavelength-dependent expression 
dMawet et al.|2005a| ): 



Mheo(^) = 



/coro(^) ^ [1 - V^]' + ^W' 

/off(^) " 



(4) 



where Icow (resp. /off) is the signal intensity when the input beam 
is centred on (resp. off) the coronagraph, €(A) is the phase error 
with respect to n, and q(A) is the ffux ratio between the two polar- 
isation components, TE and TM. Note that in the theoretical case 
where the coronagraphic profile remains identical to the original 
PSF profile, this formula of the null depth is also valid for the 
peak-to-pea k attenuation, w hich was demonstrated in an ideal 
case by Ma wet et al. ( 2005a ). 

Eq. [4] involves all the geometrical parameters of the grating 
(period A, filling factor F, depth h, sidewall angle a), illustrated 
in Fig. [2] The period A is the only fixed parameter in our design. 
It is determined by the sub wavelength limit (A/n). Considering 
a little margin for the spectral band lower bound (^ 3.4 yum), and 
considering the refractive index of diamond (^ 2.38), the period 




Fig. 2. Schematic diagram of a trapezoidal grating. The grating 
sidewalls have an angle a and an average width Fequiv^- 




FJiJing Fador F 



Fig. 3. RCWA multi-parametric simulation: mean null depth 
(logarithmic scale) over the whole L band (3.5-4. lyum) with a 
ranging from 2.7° to 3.2°. The period is set to A = 1.42 yum (SG 
limit). 



should be < 1.42yum. The other parameters of the manufactured 
grating are free parameters whose values are provided by numer- 
ical optimisation. Some of the parameters, like the depth h, the 
filling factor F and the sidewall angle a, may vary slightly during 
the etch process. The only way to measure them precisely is by 
cracking the sample along a diameter. By doing so, we measured 
that the angle a is usually 2.95° ±0.25° . Due to the uncertainty on 
the slope, we sought a robust design which performs well even 
slightly outside of the optimal parameter set. In order to find the 
target design parameters, we computed two-dimensional maps 
of the theoretical null depth as a function of the filling factor F 
and of the depth h, for several values of the angle a ranging from 
2.7° to 3.2°. The mean of all these maps (see Fig. [3]) gives the 
targeted specifications, i.e., the optimal parameters set region: 
period ^ 1.42yum, filling factor ^ 0.45, depth ^ 5.2yum, and 
sidewall angle ^ 2.95°. For these optimal parameters, the mean 
null depth over the whole L band (3.5-4.1 yum) equals 5 x 10"^. 
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Table 1. Measured profile parameters and simulated raw null 
depth for the four manufactured L-band AGPMs, and compari- 
son with the optimal design. 



Fig. 4. SEM picture of the cleaved subwavelength grating, from 
which are deduced the geometric parameters of the AGPM-L4 
profile: line width ^0.58 yum (i.e., filling factor ^0.41) and depth 
^4.7 fim. 



3.2. Manufacturing and metrology 

The AGPM pattern is etched onto a commercial optical-grade 
polycrystalline diamond window of 10 mm diameter and 300 yum 
thick, grown by chemical vapour deposition (Element Six Ltd.). 
The etching of the grooves involves techniques inherited from 
the micro-electronics industry. Electron-beam lithography was 
used to produce the pattern on a silicon wafer. Metal masking 
layers were deposited on the diamond by sputtering, and nano- 
imprint lithography was used to transfer the pattern to the top of 
the mask. Finally, reactive ion etching was used to etch the pat- 
tern first into the metal mask and then into the diamond. A sim- 



ilar manufacturing process has bee n described before ( [Karlsson 
et al.||2010 Delacroix et al.|2010 ). This process was painstak- 
ingly optimised to achieve the high pattern homogeneity, preci- 
sion and aspec t ratio necessary for half - wave plates and AGPMs 
in the N band ( [Delacroix et al.|2012a| ). This included introduc- 
ing a soft silicone stamp in the nano-imprint process and adding 
masking layers to better control the mask thickness and profile. 
For the L-band AGPM only minor adjustments had to be 
made to take into account the difficulties of etching a finer grat- 
ing. In particular, the etch rate and sidewall angle were both 
more sensitive to variations in groove width and the depth of the 
grooves was even more difficult to measure with sufficient pre- 
cision. In order to make precise measurements of the profile and 
depth, the sample had to be cracked and the cross-section im- 
aged by scanning electron microscopy (SEM), as seen in Fig.|4] 
In the N-band case, the depth could be estimated with reasonable 
precision without cracking the sample or removing the mask, 
so etching could be carried out in steps with depth measure- 
ments in between until the correct depth was reached. For the 
L-band, this was not possible and a slightly diff'erent approach 
had to be taken. After etching the metal mask, the sample was 
imaged by SEM in order to measure the line width. This width 
did not change noticeably during diamond etching, so the top 
width (FA in Fig. |2]) of the final structure was known before- 
hand. Test samples with identical masks as the AGPM sample 
were then etched and cracked in order to determine the etch rate 
and sidewall angle. Two samples etched to diff'erent depths close 
to the final depth were required, as the etch rate slows when the 
grooves get deeper. The optimal depth was then recalculated by 



name 


period A 


fillF 


depth h 


angle a 


Mheo 


AGPM-Ll 
AGPM-L2 
AGPM-L3 
AGPM-L4 


1.42 yum 
1.42 yL/m 
1.42 yum 
1.42 yL/m 


0.35 
0.36 
0.44 
0.41 


4.2 yum 
3.6yL/m 
5.8 yum 
4.7yL/m 


2.40° 
2.65° 

3.25° 
3.10° 


0.0267 
0.0116 
0.0048 
0.0010 


optimum 


1.42 yum 


0.45 


5.2 yum 


2.95° 


0.0005 




Fig. 5. SEM picture of the antireflective grating (ARG) etched 
on the backside of the AGPM. 



RCWA with the measured top width and sidewall angle fixed. 
Finally, the AGPM was etched for an appropriate length of time 
to achieve the desired depth. A detailed account of the masking 
and etching process used both for the N band and the L band will 
appear elsewhere (Forsberg et al., 2013, submitted to Diamond 
and Related Materials). The grating depth of the finished compo- 
nents was confirmed by careful parallax measurements between 
SEM pictures at five diff'erent tilt angles. 

Four finished components have been made to date. The first 
two, AGPM-Ll and AGPM-L2, were initial tests. These were 
etched with the exact same methods as N-band components. 
Further, the e-beam-written master used in the patterning of 
these two had too narrow lines. AGPM-L3 was an experiment 
with a diff'erent etch recipe giving a grating profile of trapezoidal 
walls with a triangular to ps, d esigned to reduce the problem of 
ghost reflections (see Sect. |3.4| ). AGPM-L4 finally, was produced 
using the method described above to achieve a close to optimal 
depth for its actual fill factor and sidewall angle. That we were 
able to produce nearly optimal grating parameters with so few 
trials directly reflects the maturity of the fabrication process. The 
grating parameters of all four components and their simulated 
raw null depth can be found in Table [T] (with AGPM-L3 profile 
approximated as trapezoidal). In the rest of this paper, only the 
best component, AGPM-L4, will be discussed. 



3.3. Antireflective grating 

In order to avoid incoherent reflections on the internal sides of 
the AGPM, which would not be cancelled by the optical vortex, 
an anti-reflective grating (ARG) was etched on the backside of 
the component, using a very similar diamond etching technique 
( [Karlsson & Nikolajeff|2003| with a binary square shaped struc- 
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Fig. 6. L-band transmission measurements of AGPM-L4 with an 
antireflective grating (ARG) etched on its backside. The figure 



Fig. 7. L-band expected performance of AGPM-L4, calculated 
as the sum of the theoretical raw null depth A^theo and the ghost 



also shows the calculated transmission curves, with/without the contribution A^ghost- The figure also shows the theoretical raw null 
substmte"' '^'*^'^'*°"^ *^ ^^^' ^"""^ ^°' ^ '"'°°* '^'^"'""'^ depth obtained with the optimal design presented in Table[T] 



ture as can be seen in Fig. [5] The backside reflection, which is 
high for diamond in the L band (~ 17%), is significantly reduced 
(~ 1.9%) thanks to the use of the ARG. The total transmission 
of the AGPM was measured on a Perkin-Elmer 983 spectropho- 
tometer and compared to the theoretical transmission calculated 
with RCWA (Fig.|6]). The absorption, which is significant in the 
L band (especially around 4 yum), was taken into account in the 
calculation. The result is very satisfying, with an average total 
transmission -87%. 



3.4. Expected performance 

The performance of the coronagraph, i.e., the measured raw null 
depth A^AGPM, can be described as the sum of two terms: 



A^AGPM = Mheo + A^ghost • 



(5) 



The first term is the theoretical null depth (Eq. [4]), which is lim- 
ited by the manufacturing accuracy. The second term is induced 
by the ghost, i .e., t he unwanted internal double reflection de- 
scribed in Sect.|3.3l The ghost contribution A^ghost is the ratio be- 



tween the ghost intensity and the signal intensity when the input 
beam is centred ofl' the coronagraph: 



M 



ghost 



hff TsqT ARG 



^ARG^SG 



(6) 



where Tsg, T^rg^ ^sg and Rarg are the transmit- 
tances and reflectances for both the subwavelength grating 
(SG) and anti-reflective grating (ARG) interfaces. In our case 
^ARG^SG - 10"^, which is quite small but not negligible in the 
case of the best component (AGPM-L4), whose theoretical null 
depth is also 10"^ (see Table IJ). The expected performance of 
AGPM-L4 was calculated and is presented in Fig. |7] The mean 
null depth over the whole L band (3.5-4.1 yum) is 2 x 10"^. 



4. Laboratory results 

After the manufacturing of difl'erent AGPMs with slightly vary- 
ing grating heights and filling factors, coronagraphic tests were 
carried out in the L band at the Observatoire de Meudon, on the 
same optical bench (YACADIRE) as the one used for character- 
ising the VLT-SPHERE coronagraphs ( [Boccaletti et al.|2008] ). 



4.1. Description of the coronagraphic test bench 

The setup of the coronagraphic test bench, taking into account 
some minor modifications, is represented in Fig.[8| The bench is 
approximately 1.5 m long and located in a clean room. A tung- 
sten lamp is used to feed a single-mode fibre (our star simulator) 
transmitting wavelengths up to 5 jim. Gold coated parabolic and 
flat mirrors are placed successively in the light path between the 
fibre and the detector, so that the beam alternates from pupil to 
focal planes. This choice of a complete set of reflective optical 
components avoids inherent chromatic aberrations. The AGPM 
to be tested is placed in an intermediate focal plane between two 
parabolic mirrors. With a focal distance of F = 750 mm for the 
first parabola and a beam diameter of 18.75 mm, the f-number 
of the beam impinging on the AGPM is F/40. The core of the 
Airy spot created on the AGPM is FA/D ^ 150yum wide. A Lyot 
stop blocking the diff'racted light is then located after the second 
parabolic mirror, reducing the pupil to 80% of its diameter. The 
beam is subsequently focused at F/94 with the third parabolic 
mirror on a HgCdTe infrared detector placed in a cryostat cooled 
down to 60 K. Inside the cryostat, the beam first passes through 
a wide L-band filter (3.5-4.0 yum), whose transmission character- 
istics are given in Fig. [9] 

4.2. Working and measurement principles 

We started our coronagraphic tests by placing the AGPM pre- 
cisely at the focus of the beam. A coarse positioning was 
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Tungsten 
lamp 



Shutter 




Parabolic 

mirror #3 

F = 750 mm 



Cryostat with 
HgCdTe detector 



Fig. 8. Optical layout of the YACADIRE coronagraphic test 
bench. 




3.5 3.75 

Wavelength (|Lim) 



4.25 



Fig. 9. Transmission curve measured at cold for the wide L-band 
filter (Spectrogon BBP-3500-4000nm). 



achieved using a visible laser, by looking at the diff'raction rings 
created upon propagation through the AGPM grating. The fine 
positioning was then performed in the L band, by minimising 
the flux reaching the camera. A second figure of merit in this 
optimisation was the shape of the coronagraphic PSF, which is 
expected to be symmetric when the AGPM is properly centred in 
X and y (the horizontal and vertical axes). The adjustment of the 
focus (z axis) was then based solely on the quality of the on-axis 
light extinction. 

Series of 100 frames of 40 ms were recorded with the AGPM 
placed at the optimal position (referred hereafter as the corona- 
graphic frames) and subsequently with the AGPM far from this 
position (typically at 1 mm) to measure a reference PSF with- 
out coronagraphic eff'ect, but with the beam still propagating 
through the diamond substrate (referred hereafter as the off-axis 



frames). These two sorts of frames (coronagraphic and off'-axis) 
were used to evaluate the null depth as described hereafter (Sect. 
|4.3| ). Series of 100 frames were also recorded for various off'sets 
of the AGPM from the optimal position, both in x, y and z. In 
each case, 10 to 20 background frames were recorded to allow a 
basic treatment of our images. All exposure times were limited 
to 40 ms due to the saturation of the thermal background. 

4.3. Data processing 

The first step consisted in computing the median of the back- 
ground images, before subtracting it from either the corona- 
graphic or the off'-axis frames. As we noticed that the back- 
ground was sometimes slightly over- subtracted (resulting in neg- 
ative values in the reduced images), a scaling of the background 
was implemented by using the ratio of the mean intensity in 
identical rings in both the coronagraphic (or off-axis) and the 
background frames. A flat field was measured and subsequently 
divided from the background- subtracted images. Noting /scai the 
aforementioned scaling factor, the reduced images /red j are given 
by: 



A-edj' 



^j /scaHS,med 



If 



with /seal 



<ring(/eoro)) 
<ring(/fi,med)) 



(7) 



where Ij is the j* image within the coronagraphic or off'-axis 
frames, /fi,med the median of the background images, ?> the nor- 
malised flat field, and () represents the average over a number of 
pixels. Next, bad pixels and cosmic rays were adequately elim- 
inated. Finally, all the reduced images were stacked to increase 
the signal to noise ratio of both the off'-axis and coronagraphic 
images. Aperture photometry was then carried out on these final 
images in order to measure the flux in the coronagraphic and ofl'- 
axis PSF, and thereby estimate the performance of our AGPM 
across the full L band. 

Although the peak-to-peak attenuation is a widely used met- 
rics to quantify the coronagraphic performance, the presence of 
the AGPM slightly changes the PSF profile near the axis. This 
phenomenon is characteristic of all coronagraphs. Therefore, a 
more robust metric needs to be used, taking into account both the 
coronagraphic and the off'-axis frames mentioned in Sect. |4.2[ In 
our case, the figure of merit is the raw null depth, defined as the 
ratio between the integrated flux over a certain area around the 
center of the final coronagraphic image and the integrated flux 
over the same area in the off-axis image. The question of the rel- 
evant size of this area is interesting. Using the full images would 
lead to a pessimistic result as it would integrate a lot of back- 
ground and high frequency artifacts. Using the sole central peak 
seems to be more appropriate. As we wish to com pare our re- 
sults with the theoretical predictions of Riaud & Hanot (2010]), 
we decided to use an area of radius equal to the FWHM of the 
PSF (~ A/D). This generally results in a good agreement with 
the theoretical formula (Eq.[5]), because most of the intensity in 
the image is contained within the FWHM. Our raw null depth is 
thus defined as: 



A^. 



AGPM 



Jo Jo Icow(r,0)rdrd6 



f~f:Ur,e)r,r,e 



(8) 



where /off and /eoro are the medians of the reduced off-axis and 
coronagraphic images, respectively. This definition would be 
perfectly equivalent to Eq. |5] if the coronagraphic profile was 
identical to the original PSF profile. 
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Fig. 11. Theoretical and measured off-axis transmission as a 
function of the angular separation between the beam and AGPM 
centre, expressed in resolution elements (A/D). 



Fig. 10. Azimuthally averaged PSF (plain blue curve) and coro- 
nagraphic profiles (dash green curve) with confidence interval 
(green shaded region). The top-right inset shows the off'- axis 
(left) and coronagraphic (right) images on a common linear 
scale. The off'-axis image is thresholded and the coronagraphic 
image underexposed, conveying the true raw dynamic range 
achieved by the coronagraph, with a measured raw null depth 
of 2 X 10"^ and a raw contrast of 6 x 10"^ measured at 2A/D. 



4.4. Measured performance and comparison with theory 

The raw null depth measured by the described procedure was 
2 X 10"^ over the L band (or alternatively, a rejection ratio 
R = l/N ^ 500), which is perfectly in line with the expected 
performance based on our RCWA modeling (see Sect. 



3.4). In 



Fig. [TOj the measured radial profile for the coronagraphic PSF 
(dotted line) is compared to the off'- axis PSF profile (solid line). 
These profiles were obtained by computing the azimuthal me- 
dians in concentric rings of various sizes in our images, with a 
subsequent normalisation to the peak value of the off-axis PSF. 
We computed a confidence interval on our median coronagraphic 
profile by measuring a series of 45 coronagraphic profiles on an 
integration time of about 40 sec each (i.e., 1000 frames), and by 
computing the standard deviation of the profiles in each individ- 
ual ring. The main source of ffuctuations in the coronagraphic 
profile level is associated to background instabilities. The green 
shaded region in Fig. 10 shows that the coronagraphic profile 
is measured with a sufficient signal-to-noise ratio up to about 
3A/D. Beyond this point, the intensity in the median corona- 
graphic profile is dominated by background subtraction residu- 
als. The raw contrast delivered by the AGPM at 2A/D equals 
6 X 10"^, which corresponds to 10.5 mag. Beyond 3A/D, the 
contrast quickly falls below 10"^ (i.e., 12.5 mag) and cannot be 
measured any more. 

In a second phase, we explored the dependency of the coro- 
nagraphic attenuation as a function of the offset of the AGPM 
with respect to the optimal position. The results of these mea- 
surements are shown in Fig. [TT] The agreement with the theo- 
retical off-axis transmission profile of Riaud & Hanot ( |2010j ) is 
reasonable, which further suggests that our AGPMs produce the 
textbook vortex effect. Finally, we measured the sensitivity of 



the null depth to the focus. The null depth was divided by a fac- 
tor 2 for a defocus of about ±1.5 mm at F/40, translating to a 
defocus aberration of ±120nm peak-to-valley (~ 30 nm RMS). 



5. Conclusions and perspectives 

In this paper, we have presented the design, manufacturing and 
measured performance of the first vortex phase mask corona- 
graphs based on the use of sub wavelength gratings. Our Annular 
Groove Phase Masks, etched on diamond substrates, produce on- 
axis light rejection reaching a factor up to 500 across the L band 
(3.5-4.0yum), which represents the best broadband performance 
of any mid-infrared phase mask to our knowledge. Assuming 
a high-quality input wave front, the use of our L-band AGPM 
on a telescope would result in an achievable raw contrast rang- 
ing from about 10"^ (7.5 mag) at an angular separation of lA/D 
to about 6 X 10"^ (10.5 mag) at 2A/D, and 10"^ (12.5 mag) at 
3A/D. This underlines the great potential of AGPM-L for exo- 
planet detection and characterisation at small angular separation 
from bright nearby stars. 

In November 2012, AGPM-L3 was installed at the focus 
of N AQS-CQNICA, the infrared adaptive-optics camera of the 
VLT ( [Lenzen et aTl[2003l [Rousset et aL][2Q03] ). Although the 
NACO image quality is not perfect in the L band, with a Strehl 
ratio around 70-85%, the very first on-sky tests suggest that the 
AGPM coronagraphic mode lives up to our expectations . These 
recent results will be presented in a forthcoming paper ( Mawet 
[eraLl|20T3] ). The installatio n of AGPM-L4 on LM IRCam, the 
L-band camera of the LBTI ( [Skrutskie et al.|2010i ), is currently 
on-going. Thanks to the very high Strehl ratio delivered by the 
LBT adaptive optics at L band (> 95%, [Esposito et al.|[2gTT] ), 
we expect that AGPM-L4 will unleash its full potential on this 
instrument. 

There are currently two main limitations to the performance 
of our L-band AGPMs. The first one is the presence of a ghost, 
reffected back and forth between the parallel faces of the dia- 
mond substrate. To reach starlight extinction better than 1000, 
a more aggressive antireffective solution will need to be used, 
based e.g. on an improved subwavelength structure or on a com- 
bination of the current design with an appropriate coating. The 
second limitation resides in the control of the grating param- 
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eters during manufacturing, including the slope of the grating 
walls. RCWA simulations show that, if the grating parameters 
can be made perfect, a theoretical on-axis attenuation of almost 
10""^ could be reached. It is however expected that for such deep 
extinction, local phenomena not taken into account by RCWA 
simulations could become dominant. Such considerations will 
be fully investigated when designing and manufacturing second- 
generation AGPMs, for applications on future extremely large 
telescopes. 
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